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ABSTIUCT 


This  report  presents  a  study  of  some  of  the  problems  that  are  of 
importance  in  the  design  of  a  wire-grid  lens  of  the  circularly  symmetric 
type. 

The  first  problem  investigated  is  that  of  a  double-wire  square-mesh 
grid  which  will  very  often  be  used  near  the  rim  of  the  lens  to  keep  the 
grid-to-grid  spacing  suitably  small*  Design  formulas  for  the  quasi-statie 
equivalent  dielectric  constant  of  such  grids  have  been  derived*  The 
validity  of  these  formulas  was  tested  experimentally* 

Another  problem  that  was  investigated  is  that  of  synthesizing  the 
index  of  refraction  in  a  circularly  symmetric  lens  when  the  index  of 
refraction  is  specified  in  part  of  the  lens  and  when  it  is  required  that 
the  rays  leaving  the  lens  must  form  a  collimated  beam  in  as  large  an 
aperture  as  possible*  Design  formulas  representing  the  solution  of  this 
problem  have  been  derived*  These  formulas  are  most  useful  for  wire-grid 
lens  designs  because  the  grid-to-grid  spacing  in  the  outer  part  of  the 
lens  is  often  determined  by  mechanical  considerations  rather  than  by 
electrical  considerations* 
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1. 


INTRODUCTION 


In  recent  years »  there  has  been  a  need  for  scanning  directional 
antennas  which  can  be  used  in  the  HF  frequency  band  (3-30  Me)  for 
direction-finding  applications^  Artificial  dielectrics  practical  and  economical 
at  these  frequencies  have  not  existed,  however.  A  most  important  structure, 
which  makes  possible  lens  designs  in  the  HF  band,  is  a  pair  of  wire  grids 
with  either  square,  hexagonal  or  triangular  meshes.  A  pair  of  wire  grids 
can  be  used  for  lens  designs  because  the  velocity  of  a  wave  propagating 
between  the  grids  can  be  controlled  by  changing  the  ratio  between  the  grid 
spacing  and  the  mesh  size.  The  wire-grid  lens  as  a  direction-finding  antenna 
and  problems  related  to  the  design  of  such  antennas  is  the  subject  of  investi¬ 
gation  in  this  report. 

A  great  number  of  problems  have  to  be  solved  to  realize  a  wire-grid 

lens.  Part  of  these  problems  have  already  been  solved  previously  through 

support  by  the  United  States  Air  Force  under  Contracts  AF  19(604)-2240  and 
12  3 

AF  19(b04)-8059.  ’  '  The  work  on  these  contracts  was  carried  out  at 
Stanford  Research  Institute  and  included  a  determination  of  the  equivalent 
static  dielectric  constant  of  a  pair  of  grids  with  square,  hexagonal  and 
triangular  mesh  in  relation  to  the  grid  dimensions.  Also  investigated  were 
the  dynamic  properties  of  a  pair  of  grids  by  an  experimental  method.  The 
usefulness  of  the  results  obtained  was  demonstrated  by  testing  an  experimental 
model  of  a  wire-grid  lens  designed  according  to  the  developed  design  informa¬ 
tion. 
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While  the  investigations  of  the  wire -grid  lens  undertaken  at  Stanford 
Research  Institute  led  to  important  basic  data  needed  in  a  practical  design^ 
many  problems  yet  have  to  be  solved  before  a  practical  wire -grid  lens  can 
be  realized.  Among  these  problems  are: 

(1)  An  extended  theoretical  investigation  of  the  dispersion 
and  anisotropy  of  a  pair  of  parallel  wire  grids 

(2)  An  investigation  of  grid  structures  whose  equivalent 
dielectric  constant  can  be  controlled  continuously, 

not  only  by  changing  the  grid  spacing  but  also  be  changing 
the  mesh  structure;  for  example,  by  changing  the  effective 
mesh  size 

(3)  An  investigation  of  various  radiating  structures  to 
provide  the  required  matching  between  the  lens  and 
free  space 

(4)  An  investigation  of  various  feed  structures 

(5)  An  investigation  of  circularly  symmetric  lenses  whose 
dielectric  constant  can  be  specified  in  part  of  the  lens 

((^  An  investigation  of  the  radiation  pattern  of  a  wire -grid 
lens  including  the  effect  of  dispersion  and  anisotropy. 

In  this  report,  problems  (2)  and  (5)  are  treated. 

In  one  of  the  reports^  referred  to  above,  square -mesh,  hexagonal - 
mesh  and  triangular -mesh  grids  were  investigated.  As  a  result  of  the  investiga- 
tion,  it  wa«  found  that,  for  the  same  mesh  area,  hexagonal -me  eh  gride  exhibit 
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less  anisotropy  and  dispersion  than  both  square-mesh  grids  and  triangular- 
mesh  grids.  From  an  economical  point  of  view,  however,  the  square  mesh 
will  most  often  be  superior  to  the  hexagonal  mesh  because  a  square-mesh 
grid  is  much  simpler  to  construct  than  a  hexagonal-mesh  grid.  In  this 
report,  therefore,  only  square-mesh  grids  are  investigated. 
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2. 


DOUBLE-WIRE  SQUARE-MESH  GRIDS 


The  equivalent  dielectric  constant  of  a  pair  of  square-mesh  grids 
is  a  function  of  the  ratio  between  the  grid  spacing  and  the  mesh  size.  The 
smaller  this  ratio^  the  larger  the  dielectric  constants  For  very  closely 
spaced  grids,  the  dielectric  constant  approaches  2  while  it  approaches  1 
for  grids  that  are  far  apart  relative  to  the  mesh  size.  At  the  rim  of  a 
circularly  symmetric  wire-grid  lens,  it  is  usually  necessary  that  the 
equivalent  dielectric  constant  be  close  to  the  value  L  This  requires  the 
spacing  between  the  grids  to  be  very  large  near  the  rim  of  the  lens.  In 
many  cases,  the  required  spacing  between  the  rim  of  the  lens  turns  out 
to  be  larger  than  practically  desirable.  To  illustrate  this,  assume  that 
a  dielectric  constant  of  1,  1  is  desired  at  the  rirn  of  the  lens.  For  a 
square -mesh  grid  lens  with  a  constant  mesh  size  of  5  feet,  which  would 
be  a  reasonable  mesh  size  for  an  HF  lens,  it  is  then  necessary  to  use  a 
spacing  of  about  35  feet.  This  spacing  could  be  reduced  by  a  factor  of 
between  Z  and  3  by  reducing  the  mesh  size  by  a  factor  of  Z.  However, 
reducing  the  mesh  size  by  a  factor  of  Z  increases  the  cost  of  the  lens 
by  nearly  a  factor  of  Z  and  is  therefore  not  a  practical  solution,  A  practical 
solution  to  the  problem  was  proposed  by  Dr,  Ro  L*  Tanner.  TRG-West, 
prior  to  the  contract  period;  by  using  a  double-wire  square -mesh  grid, 

(see  Fig,  1)  with  varying  spacing  d  between  the  wires  of  the  double-wire 
near  the  rim  of  the  lens,  the  effective  mesh  size  can  be  reduced  by  a  factor 
of  Z  near  the  rim  of  the  lens  while  keeping  the  mesh  size  unchanged  in  the 
larger  part  of  the  lens. 


-4- 


In  this  sectioHf  formulas  will  be  derived  for  the  equivalent  dielectric 
constant  of  a  pair  of  double -wire  square -mesh  grids  in  dependence  of  the 
ratio  between  the  grid  spacing  and  the  mesh  size,  and  of  the  wire  radius 
and  of  the  spacing  between  the  wires  of  the  double  wire.  This  represents 
the  basic  information  needed  for  designing  a  pair  of  double -wire  square -mesh 
grids«  The  dielectric  constant  will  be  derived  by  a  method  analogous  to  the 
method  which  was  previously  used  to  evaluate  the  dielectric  constant  of  a 
pair  of  single -wire  square -mesh  grids.^ 

The  equivalent  static  dielectric  constant  can  be  expressed  as  the 
ratio  between  the  capacitance  per  cell  of  the  grids,  and  the  capacitance  per 
cell  of  the  grids  with  one  set  of  parallel  wires  removed*  The  capacitance 
will  be  calculated  from  the  variational  expression 


where  p  is  the  charge  stored  per  unit  length  of  the  wires,  and  (p  is  the 
electrostatic  potential  produced  at  the  wire  surface.  When  p  is  the  correct 
charge  distribution,  p  is  constant  along  the  wires.  Otherwise,  p  changes  on 
the  wire  surface.  However,  since  Eq,  (1)  is  stationary  with  respect  to  small 
changes  of  p  ,  this  equation  can  provide  a  reasonably  accurate  expression 
for  C  even  if  p  is  not  too  close  to  the  correct  distribiition.  Since  the  capaci¬ 
tance  can  be  easily  evaluated  for  a  constant  charge  density  on  the  wires,  and 
this  charge  distribution  should  be  reasonably  accurate  except  near  the  points 
where  wires  intersect  each  other  (where  the  charge  density  drops  to  a  lower 
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value)  we  shall  now  assume  that  the  charge  is  uniformly  distributed  on 


the  wires. 

For  a  constant  charge  density  we  can  take  advantage  of  the  fact 
that  the  electrostatic  field  in  a  cell  is  a  sum  of  two  dimensional  fields* 
Accordingly,  we  can  use  the  theory  of  complex  functions  to  find  the  potential 
of  the  line  charges  in  the  cell, 

A  single  cell  of  the  grids  is  shown  in  Fig.  2,  A  cross  section  of 
this  cell  is  shown  in  Fig.  3  y  =  0  and  y  =  b  are  magnetic  walls  at  which 
the  normal  derivative  of  the  electrostatic  potential  has  to  be  zero.  Consider 
first  the  potential  produced  by  the  line  charges  parallel  to  the  z-axis.  Assuming 
the  xy-plane  to  be  a  complex  z -plane  (z  =  x  +  iy),  transform  the  strip  0  <  y  <  b 
to  the  upper  half  of  a  -plane  by  the  transformation 


(2) 


The  magnetic  walls  y  =  0  and  y  -  b  thereby  transform  to  the  real  axis  of  the 
{^-plane.  The  positions  of  the  line  charges  transform  as  follows 
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Fig.  3  Transformation  of  Grid  Cell  in  Complex  Plaiu* 
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The  positions  are  marked  in  Fig.  3.  Imaging  the  four  line  charges  in  the 
upper  half  of  the  (^-plane  in  the  4 -axis,  the  system  of  charges  shown  in 
Fig.  3B  is  obtained. 

Assuming  that  the  charge  per  unit  length  of  the  wires  is  q,  the 
complex  potential  in  the  4~plAne  is  found  to  be 


I 


L(r-CXfH.vl 
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where  c  is  the  dielectric  constant  of  the  medium, 
o 

The  real  potential  in  the  z -plane  is  the  real  part  of  4.  Intro¬ 
ducing  Eq.  (?.)  above,  we  find  then 


(8) 
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The  potential  at  the  plane  x  =  0  mid  between  the  grids  (this  is  an  electric 
plane  because  of  symmetry)  is  found  from  this  expression  to  be 


a,  .  ^ 

f  (ground  ^ 


(9) 


We  want  the  potential  at  x  =  0  to  be  zero.  This  can  bo  achieved  by  subtracting 
£q.  (9)  from  Eq.  (8).  Doing  so.  we  find 


(10) 


To  evaluate  the  capacitance  according  to  Eq.  (1),  we  need  to  know 

the  electrostatic  potential  at  the  wire  surface,  this  is,  at  a  distance  r  from 

o 

the  points  or  According  to  Eq.  (lO),  this  potential  is  very  nearly 

constant  and  equals 


=  S- 
Xire, 


-{in 


[  l-06\ 

iiv)  1  1 

VCoshijW] 

1 

(11) 


-11- 


The  total  charge  on  two  parallel  wires  in  a  cell  of  the  double -wire 
grids  is  2  qb.  The  capacitance  of  one  set  of  parallel  wires,  introducing 
Eq«  (11)  in  Eq*  (I),  therefore  is 


(12) 


The  capacitance  of  one  grid  cell  with  both  sets  of  parallel  wires  in 
place  can  be  expressed  as 


(13>y 


where  the  integration  is  extended  over  one -half  the  length  of  one  of  the  wires 
carrying  the  charge  q  per  unit  length.  The  potential  4>  equals  the  sum  of 
the  above  constant  potential  (see  Eq.  (11))  and  a  varying  potential  4)^ 
which  is  produced  by  the  charge  on  the  set  of  wires  perpendicular  to  the 
set  of  wires  producing  the  constant  potential  Equation  (13)  reduces  to 


c  = 
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<b  is  evaluated  along  the  line  z  =  a  -  r  +  iy  (sec  Fig.  3).  Introducing  this 
value  of  z  in  Eq.  (8).  wc  find 


(15) 


This  can  be  reduced  to 


[cc6h  ■'■  Cos 
-/n[cosi 

-/&t  /  Cesh  (r^ ^]j 
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When  thie  expreeaion  it  integrated  over  y  from  y  s  0  to  y  «  >  ,  we  find 

S<p,dH'-  f(a~r)  (17) 

c 

Introducing  this  in  Eq.  (14),  the  following  expression  is  found  for  the  capaci¬ 
tance  of  a  single  cell  of  a  double -wire  grid. 


The  equivalent  static  dielectric  constant  of  a  pair  of  double -wire 
square -mesh  grids  is  therefore  given  by  the  following  expression. 


t 


(19) 


where  G  is  defined  by 
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(20) 


I  I  ~  CdsI^^SJ!) 


«  i.  .a.i.,  .„„  E,.  „„  ««.  .  -  3 

1  for  courre,  expected. 

Tte  eccorec,  of 

. .  °'  *  »"  <He  ^re.  „e,  „o.  Oe  ver, 

I  od  cloee  to  the  point.  »here  the  grid  »ire.  intereect.  Yet  it  ™ay  be 

improve  the  accuracy  of  E,.  (19)  by  conaidering  a  charge 

dl.tribution  that  approximate,  the  correct  chare  ,  .  c 

distribution  bettor  than 

the  conetant  charge  den.ity.  The  math.maticallv  .i  . 

n^^“caUy  simplest  charge  dis¬ 
tribution  which  will  give  an  ir^ 
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rr::; "  * 

*  found  to  be  expressed  by 
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(21) 


’  ha,<^a/ 


where  c  is  the  dielectric  constant  defined  in  Eq.  (19/i  and  cC  is  related  to  the 
functions  a.  and  a  by 

A  M 


o<r- 


a~Z 

ZQ-a, 


(22) 


while  a^  and  a^  are  related  to  the  grid  dimensions  and  to  the  function  G  (Eq.  20) 
as  follows 


(23) 


and 


1 

^  G*¥Jrf 


(24) 
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Equation  (21)  has  been  plotted  in  Fig.  4  which  shows  c  as  a 

function  of  r  ior  various  values  of  ^  ,  and  for  ~  =  200  which  value  is 
b  b  r 

o 

close  to  the  value  used  in  the  experiments  described  below. 

To  test  the  validity  of  Eq,  (21),  experimental  measurements  were 

made  using  a  rectangular  resonator  operating  in  the  TE^^^  modes  in  which 

the  top  wall  of  the  resonator  was  replaced  by  the  wire  grid.  The  method 

of  determining  the  grid  dielectric  constant  from  the  resonant  frequencies 

of  tha^ resonator  was  described  in  a  previous  report,  ^  The  agreement 

between  theory  and  experiment  was  excellent.  For  a  grid  resonator  with 

the  parameters  ^  =  0,251,  0,124  and  — =  230,  a  value  of  1.  4Gl»  was 

b  D  r 

o 

measured  for  the  equivalent  dielectric  constant.  The  theoretical  value 

predicted  by  Fig,  4  is  1,  395,  For  another  grid  resonator  described  by 
ad  b 

the  parameters  —  =  0.  6,  —  =  0.  124  and  —  =  230,  a  value  of  1,  21  was 
DO  r 

o 

measured  for  c.  The  theoretical  value  is  1.  213  according  to  Fig.  4,  It 
can,  therefore,  be  concluded  that  the  accuracy  of  Eq,  (21)  is  very  good. 
Equation  (21)  was  programmed  as  a  sub-routine  for  an  electronic  computer. 
This  sub-routine  will  baused  in  subsequent  computations  on  wire«-grid  lenses. 
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3.  CIRCULARLY  SYMMETRIC  LENS  WITH  PARTIALLY 

PRESCRIBED  INDEX 

3.  1  General 

For  many  applications  of  a  circularly  symmetric  lens  using  an 
artificial  dielectric  for  the  lens  material,  the  variation  of  refractive  index 
required  for  a  Luneburg  Lens  is  not  practical.  In  particular,  a  maximum 
index  of  refraction  less  than  'y/~2  is  desirable  for  a  wire -grid  lens.  One 
method  by  which  the  maximum  index  of  refraction  required  to  collimate 
the  rays  leaving  the  front  part  of  the  lens  (opposite  to  the  feed)  can  be 
made  less  than  "V^ is  to  place  the  feed  a  certain  distance  away  from  the 
rim  of  the  lens.  The  lens  with  the  displaced  feed  can  be  considered  a  rim- 
fed  lens  whose  index  of  refraction  has  been  assigned  the  value  1  in  an  outer 
ring. 

In  practice,  the  index  of  refraction  in  the  outer  part  of  a  circularly 
symmetric  lens  will  often  be  determined  by  mechanical  considerations  more 
than  by  electrical  considerations,  and  it  may  not  be  possible  to  produce  any 
prescribed  variation  of  the  index  of  refraction.  In  view  of  this  and  in  view 
of  the  above  considerations  concerning  the  maximum  index  of  refraction,  a 
practical  design  requires  the  solution  of  the  general  problem  of  synthesizing 
a  lens  whose  index  of  refraction  is  specified  in  an  outer  ring  and  which  is  fed 
by  a  point  source  somewhere  in  this  ring.  The  purpose  of  the  investigation 
in  this  section  is  to  pursue  this  problem. 
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3,  2  The  Rays  in  the  Lens 

A  cross  section  of  one -half  of  a  spherically-symmetric  lens  and 

the  path  of  a  ray  through  this  lens  is  shown  in  Fig,  5.  The  index  of  refraction 

is  specified  in  the  ring  r  <  r  <  r  .  The  point-source  feed  is  placed  at  a 

C 

distance  r,  from  the  lens  center,  r  <  r,  <  r  .  We  want  to  determine  the 
f  a  f  c 

index  of  refraction  in  the  inner  core  of  the  lens  (r  <  r  )  such  that  the  rays 
leaving  the  front  part  of  the  lens  (opposite  to  the  feed)  are  collimated  in  as 
large  an  aperture  as  possible.  To  do  so,  we  shall  first  establish  the  general 
relation  between  the  radially  varying  index  of  refraction  n(r),  the  ray  angle 
and  the  angle  of  the  turning  point,  that  is,  that  point  of  the  ray  which  is 
closest  to  the  center  of  the  lens  (see  Fig.  5). 

The  rays  in  the  lens  move  along  paths  obeying  Snell's  law  which 
can  be  written 


rtTsitjv^x  <25) 

where  ^  is  a  constant  characteristic  of  the  particular  ray.  X  is  different 
for  different  rays  and  is  related  to  the  ray  angle  4<  a*  follows 

sin\j/  (26) 

where  n^  is  the  index  of  refraction  at  the  point  where  the  feed  is  placed. 

may  also  be  related  to  the  index  of  refraction  n  and  radius  r  at  the 

o  o 

turning  point  as  follows. 
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(27) 


^  ^  r 

o  o 

The  polar  angle  of  the  turning  point  is  most  easily  found  by 
considering  the  differential  form  of  Snell's  law 


where  the  upper  sign  is  used  when  the  ray  approaches  the  turning  point  while 
the  lower  sign  is  used  when  the  ray  recedes  away  from  the  turning  point. 
When  Eq.  (28)  is  solved  with  respect  to  we  find 


(29) 


Integrating  with  respect  to  r,  we  find  the  following  expression  for  the  polar 
angle  ^ 


(30) 
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when  the  ray  is  approaching  the  turning  point.  When  the  ray  is  receding 
from  the  turning  point,  the  following  expression  should  be  used  for  ^ 


where  is  the  polar  angle  of  the  turning  point 


According  to  Eq.  (31).  the  polar  angle  4  of  the  point  where  the  ray 


leaves  the  lens  is 


The  angle  6  between  the  direction  of  the  ray  after  leaving  the  lens  and  the 
principal  direction  of  the  lens  according  to  Eq.  (25)  and  Fig.  5  is  defined 
by  the  relation 


s(k[7T’‘(§'t’9)\  *  X 


-23- 


from  which  it  follows  that 


9-^  7T -  Arcsin(^^  os) 

The  above  formulas  are  useful  for  determining  ray  paths  when  the 
.ndex  of  refraction  is  specified  throughout  the  lens.  These  formulas  form 
die  basis  for  synthesizing  a  lens  whose  index  of  refraction  is  unknown  in 
3art  of  the  lens  and  for  calculating  radiation  patterns. 
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3,  3  Synthesis  of  Lens 


We  shall  now  consider  the  synthesis  problem:  given  the  index  of 

refraction  for  r  <  r  <  r  ,  to  find  the  index  in  the  inner  region  0  <  r  <  r 
21  C 

such  that  those  rays  that  pass  through  the  inner  region  form  a  collimated 
beam  in  as  large  an  aperture  as  possible  upon  leaving  the  lens.  To  solve 
this  problem,  consider  Eq.  (32).  This  is  an  integral  equation  from  which 
the  index  n  can  be  found  when  (|>^  is  given  as  a  function  of  ^  ,  According 
to  Eq.  (26),  is  defined  in  the  range  0  <  ^  <  n^  r^.  A  certain  part  of 
this  range  belongs  to  rays  that  do  not  penetrate  the  inner  region  0  <  r  <  r 

A 

of  the  lens.  The  value  of  (]>^  {^)  in  this  part  of  the  total  range  of  ^  can 
be  calculated  from  the  given  index  of  refraction  for  r  >  r  by  Eq.  (32). 
However,  some  care  has  to  be  exercised  in  the  use  of  Eq.  (32)  if  nr  is 
not  a  monotonic  function  of  r. 

If  nr  is  a  monotonic  function  of  r  (see  Fig.  6A),  Eq.  (32)  can  be 

used  with  no  difficulties.  The  rays  that  do  not  penetrate  the  inner  region 

0  <  r  <  r  of  the  lens  are  in  that  case  those  for  which  has  a  value  in 
a 

the  range  n  r  <  ^  <  n  r  .  Physically,  the  monotonic  variation  of  nr 

&  2L  X  1 

with  r  implies  that  one,  and  always  one,  ray  passes  through  every  point 
in  the  lens. 

Assume  now  that  nr  is  not  a  monotonic  function  of  r,  as  in  the 
examples  with  two  extrema  of  the  rn*curve  shown  in  Figs.  6B  and  6C« 

In  these  figures,  r"  is  the  radius  for  which  nr  assumes  the  minimum 
while  r*  <  r"  is  the  radius  for  which  nr  has  the  same  value  as  for  r  s  r". 
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Consider  first  Fiij*  6B  where  n  r  is  the  lowest  value  of  nr  for  r  >  r  . 

°  a  a  a 

In  that  case  r  in  Eq.  (32)  cannot  be  allowed  to  have  any  value  in  the  range 
o 

from  r  to  r,,  but  the  interval  r'  <  r  <  r"  must  be  excluded  from  that  range, 
a  f  o 

Since  r^  <  r'  <  r",  the  rays  that  do  not  penetrate  the  inner  region  0  <  r  <  r^ 
of  the  lens  are  again  those  for  which  'XI,  has  a  value  in  the  range  n  r  <  ^ 

3l 

<  n^  r^.  The  physical  consequence  of  the  type  of  variation  of  nr  shown 
in  Fig,  6B  is  that  a  ray  will  not  pass  through  any  point  of  the  lens.  This 
is  illustrated  in  Fig,  7.  This  variation  is  seen  to  have  the  effect  of  producing 
regions  of  "darkness"  in  the  lens.  This  may  be  a  desirable  feature  for  many 
applications. 

In  Fig.  6C  where  r'  <  r  <  r",  the  rays  that  do  not  penetrate  the 

inner  region  0  <  r  <  r  of  the  lens  are  those  for  which  ^  has  a  value  in  the 

range  n*'r”  <  ^  <  r^.  In  contrast  to  the  nr-curves  in  Figs.  6A  and  6B, 

the  nr-curve  in  Fig.  6C  is  not  allowed  because  it  matkes  it  impossible  to 

determine  the  refraction  index  in  the  region  0  <  r  <  r  .  To  solve  Eq,  (32) 

for  the  index  of  refraction  in  the  inner  part  of  the  lens  from  a  prescribed 

variation  of  with  ^  for  the  rays  that  penetrate  the  inner  part  of  the  lens, 

assume  that  n  r  is  the  smallest  value  of  nr  for  r  >  r  ,  The  range  of 
a  a  a 

for  the  rays  penetrating  the  region  0  <  r  <  r  is  therefore 

a 

0^x.<\r^  (36) 
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Fig.  7  Ray  Paths  in  a  Lens  with  an  nr  Distribution 
as  shown  in  Fig.  bB 
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Now  separate  from  the  integral  in  Eq.  (32)  the  part  for  *  ^ 
which  is  known.  We  find  then 


(37) 


where 


(38) 


is  a  known  function  of  ^  since  ^JiOC)  will  be  required  to  be  specified  such 
that  all  rays  passing  through  the  inner  part  0  <  r  <  r  of  the  lens  form  a 

2t 

collimated  beam  upon  leaving  the  lens.  Obviously,  is  therefore  found 

from  Eq.  (35)  by  putting  6=0.  Introducing  Eq.  (33)  in  Eq.  (35),  we  find 


T 

T 

2L  [ 

f 


Let  us  introduce  this  in  Eq.  (38).  4'q(^)  1*  then  found  to  be  given  by 


(39) 


- 1 -i4rc£h(f)-  f 


dr 


(40) 
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4  5 

The  solution  of  Eq.  (37)  has  been  found  by  Kay  and  Morgan  for 

the  case  r  s  1  and  n  r  s  1.  The  solution  for  this  case  is  given  by 
a  a  a  ©  / 


7>=  f  t 


dx 


(41) 


where  4  a  parajneter  defined  by 


(42) 


These  two  equations  serve  to  determine  n  and  r. 

The  solution  of  Eq.  (37)  can  be  deduced  from  Eq.  (41)  by  trans¬ 
forming  Eq.  (37)  as  follows 


Obviously,  this  transformation  normalizes  both  r  and  nr  to  the  value  1. 

Accordingly,  the  solution  of  Eq.  (37)  can  be  deduced  from  Eq,  (41)  by 

replacing  r  by  —  ,  n  by  —  and  ^  by  ?Cln  r  ,  Substituting  also 
r  n  a  a 

a  a 

t  -  tin  r  ,  we  find 

a  a 
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'cl 


Af  i 


(44) 


and 


7} 


(45) 


The  method  used  to  derive  Eq.  (41)  restricts  the  solution  for  nr  to  be  a 
monotonic  function  of  r. 

For  the  evaluation  of  Eq.  (44),  we  have  to  solve  the  integrals  in 

Eq.  (40).  These  integrals  cannot  be  reduced  to  a  closed  form  except  in 

very  special  cases,  among  these  the  case  where  n  is  constant.  To  avoid 

having  to  integrate  numerically  first  over  r  according  to  Eq.  (40),  and  then 

over  OC  according  to  Eq.  (44),  we  shall  assume  that  the  region  r^  <  r  <  r 

of  the  lens  consists  of  a  large  number  of  layers,  each  layer  with  a  constant 

index  of  refraction.  The  resulting  stepped  index  curve  can  be  made  to 

approximate  the  actual  continuous  index  curve  as  well  as  desired  by  choosing 

a  sufficiently  large  number  of  layers.  Assume  that  the  region  r  <  r  <  r. 

consists  of  P,  layers  with  refraction  indices  n,.,  n.,,  n._,  n.  _ 

1  11  12  13  1,P, 


and  with  radii  r^^,  r^^t 


<  r 


l.Pj+l* 


-31- 


We  find  then 


Assume  similarly  that  the  region  r^  <  r  <  r^  consists  of  layers 


with  indices  n^^,  n^^,  n^^, 


p  and  with  radii 


"a,  p,ti  h,  p,ti- 


analogously 


Now  introducing  Eqs.  (46)  and  (47)  in  £q.  (40),  and  Eq.  (40)  in 
Eq.  (41),  and  introducing  also  the  integral  function 


)C 

x)-  / 


krc$in(^) 


(48) 
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which  has  been  expressed  as  a  power  series  in  Appendix  a,  we  iina 


(49) 


where 


Pz 


1 


-i  tU  0^  ,  %T,)  -1 0,  >%r^)] 


(50) 


The  numerical  evaluation  of  the  integral  in  Eq.  (48)  is  discussed 
in  Appendix  A. 

Equation  (49)  was  programmed  for  an  electronic  computer.  The 
index  of  refraction  in  this  formula  has  been  assumed  to  be  specified  in  the 

outer  oart  r  <  r  <  r  of  the  lens.  Since  the  formulas  derived  above  will  be 

a  c 
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used  to  design  a  wire -grid  lens,  instead  we  shall  specify  the  grid-to-grid 

spacing,  the  mesh  size,  and  the  spacing  between  the  double  wires  in  a  double - 

wire  square  mesh*  The  index  of  refraction  for  r  <  r  <  r  was  then  evaluated 

a  c 

by  a  sub-routine  program  for  Eq.  (21).  The  result  of  the  computations  is  a 

table  of  the  index  of  refraction  versus  r  according  to  Eq.  (49).  Use  of 

Eq.  (21)  then  gives  the  grid-to-grid  spacing  in  the  entire  lens. 

An  example  of  the  use  of  the  computer  program  for  synthesizing 

a  wire -grid  lens  shall  be  presented*  Assunie  that  the  lens  is  divided  into 

four  regions  as  illustrated  in  Fig,  8.  In  the  region  0  <  r  <  r  the  index  is 

a 

unknown.  The  mesh  size  used  in  this  region  is  b.  In  the  region  r  ^  r  <  r, 

a  b 

a  double -wire  grid  is  used  with  linearly  tapered  spacing  between  the  double 

wires.  The  double -wire  grid  provides  a  continuous  ’’match’*  between  a  single - 

wire  grid  with  mesh  size  b  and  a  single -wire  grid  with  mesh  size  In  the 

regions  r^  <  r  <  r^  and  r^  <  r  <  r^  a  single -wire  grid  with  mesh  size  ^  used. 

The  point  source  feed  is  placed  at  r  =  r^.  The  values  chosen  in  this  example 

for  r  ,  r  ,  r  and  r  ,  and  for  the  values  of  ^  and  —  at  these  radii  are  shown 
able  b  a 

in  Fig,  8.  Both  ^  and  assumed  to  vary  linearly  with  r  (except,  of  course, 

in  the  region  0  <  r  <  r  where  ^  is  unknown).  The  wire  radius  was  chosen  such 

a  b 

that  — =  1000. 
r 

o 

The  results  of  the  synthesis  for  the  lens  example  shown  in  Fig,  8  are 

plotted  in  Figs,  9  and  10,  Figure  9  shows  the  nr-product  and  Fig,  10  the 

index  of  refraction  for  the  entire  lens.  From  Fig,  10  it  will  be  observed  that 

the  index  rises  quite  sharply  near  r  s  r  =0,  75,  As  a  result  of  this,  the  nr- 

a 

curve  in  Fig.  9  exhibits  a  minimiun  for  r  >  r  .  It  was,  therefore,  necessary 

a 

to  further  specify  a  constant  index  in  a  narrow  ring  r'  <  r  <  r  such  that  the 

a 
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value  of  nr  at  r  s  r'  would  be  the  smallest  value  of  nr  for  r  >  r'.  This 
condition  was  satisfied  for  r'  =  0.  745  in  the  example  considered.  The 
computer  is  programmed  to  introduce  the  extra  ring  with  a  constant  index 
of  refraction  when  necessary. 

The  validity  of  the  computer  program  for  synthesizing  a  circvilarly 
symmetric  lens  was  checked  by  synthesizing  the  well-known  Limeburg  Lens. 
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CONCLUSIONS 


The  investigation  of  a  double -wire,  including  the  special  case  of 
single-wire,  square-mesh  grid  lens  whose  shape  and  mesh  structure  is 
specified  in  an  outer  ring  of  the  lens,  has  led  to  a  computer  program  which 
determines  the  grid  shape  in  the  entire  lens  such  that  the  rays  leaving  the 
lens  form  a  collimated  beam  in  as  large  an  aperture  as  possible. 
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APPENDIX  A 


For  •ynthesizing  a  circularly-symmetric  lens  whose  refractive 
index  is  specified  in  part  of  the  lens,  we  need  to  evaluate  the  following 
integral 


This  integral  cannot  be  expressed  by  known  functions,  and  a  numerical  evalua 
tion,  when  high  accuracy  is  required,  can  become  qmte  time  consuming  since 
the  integrand  has  a  singularity  at  the  lower  integration  limit.  Fortunately, 
however,  the  difference  1(4, a  ,x)  -  1(4,  a., x)  between  two  integrals  with 
slightly  different  values  of  a  enters  the  synthesis  problem  many  more  times 
than  the  single  integral.  The  integral  difference  can  be  evaluated  very  fast 
and  with  high  accuracy  by  expanding  the  integral  in  a  power  series  of  a. 
Expanding  from  a  =  a^  where 


a 


z 


we  find  then 
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Most  often,  terms  of  higher  order  than  3  in  ~  can  be  neglected. 

The  first-  and  third-order  derivatives  of  I  can  be  evaluated  in 
closed  form.  Differentiating  I  with  respect  to  a,  we  find 

X 

fta  “  'aJ 

*  -/trcsk^l-Z  -fr-f)] 

Differentiating  this  expression  another  two  times,  the  following  expression 
is  found  for  the  third -order  derivative 


/J{Ux) 

da? 


ZArcth 


r- 
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when  it  is  required  that  the  rays  leaving  the  lens  mast 
form  a  collimated  beam  in  as  large  an  apertnre  as  possible., 
Derived  design  formalas  representing  the  solution  of  this 
problem  are  most  osefnl  for  wire-grid  lens  designs  because 


